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Abstract 

Probabilities of crossing on same-spin clusters, seen as order parameters, have been 
introduced recently for the critical 2d Ising model by Langlands, Lewis and Saint- 
Aubin. We extend Cardy's ideas, introduced for percolation, to obtain an ordinary 
differential equation of order 6 for the horizontal crossing probability tt^. Due to the 
■^j- ■ identity 7T/,(r)+7r/ l (l/r) = 1, the function tt^ must lie in a 3-dimensional subspace. New 

measurements of irh are made for 40 values of the aspect ratio r (r £ [0.1443,6.928]). 
| These data are more precise than those obtained by Langlands et al as the 95%- 

confidence interval is brought to 4 x 10~ 4 . A 3-parameter fit using these new data 
determines the solution of the differential equation. The largest gap between this 
solution and the 40 data is smaller than 4 x 10 -4 . The probability nhv of simultaneous 
horizontal and vertical crossings is also treated. 
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Les probabilites de traversee sur les plages de spins identiques, vues comme para- 
. metres d'ordre, ont ete introduites recemment pour le modele d'Ising bidimension- 

nel critique par Langlands, Lewis et Saint-Aubin. Nous etendons les idees de Cardy, 
presentees pour la percolation, afin d'obtenir une equation differentielle ordinaire du 
sixieme ordre pour la probabilite d'une traversee horizontale ir^. II suit de l'identite 
T^hir) + 7>7i(l/ r ) = 1 doit reposer dans un sous-espace de dimension trois. De 

nouvelles mesures de iTh sont faites pour 40 valeurs du rapport r des cotes du rectangle 
(r 6 [0.1443,6.928]). L'intervalle de confiance de 95% est ici reduit a 4 x 10~ 4 . La 
solution de l'equation differentielle est choisie en minimisant l'ecart quadratique entre 
la solution et ces nouvelle donnees. Le plus grand ecart entre les donnees et la solution 
ainsi obtenue est plus petit que 4 x 10 -4 . La probabilite de traversees horizontale et 
verticale simultanees est egalement etudiee. 

Short Title: Crossing probabilities in the 2d Ising model 
PACS NUMBERS: 0.50.+q, 64.60.Fr, 64.60.Ak, 11.25.Hf 
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1 Introduction 



The probability 7r^(r) of crossing on open sites inside a rectangle of aspect ratio r has been 
measured at p c for several models of percolation in and |7]]. These simulations support 
hypotheses of universality and conformal invariance of this function 71^ (r) and of several 
others. Cardy's contemporaneous work |I| offered a prediction for 717, using conformal field 
theory. His analytic expression agrees with the simulations within statistical errors and 
provides further support that these crossing probabilities are order parameters with the 
usual critical properties. 

This paper presents similar evidence, both numerical and analytic, for crossing probabil- 
ities on same-spin clusters of the two-dimensional Ising model at criticality. These crossing 
probabilities are not traditional order parameters for the Ising model. It is not a priori obvi- 
ous that they are not identically zero or one in any dimension d. In dimension two however 
the simulations carried out in || indicates clearly that they are non-trivial functions and 
that they are likely to satisfy the same hypotheses (or even more restrictive ones) of uni- 
versality and conformal invariance. Moreover physical quantities for which conformal field 
theory gives quantitative predictions that can be readily verified by simulation are always 
welcome. 

Let a triangular lattice be oriented in such a way that sites on horizontal lines 
distance of one mesh unit. A rectangle of height V and width H is superimposed on the 
lattice. The width H is measured in mesh units but V is the number of horizontal lines in 
the rectangle. A configuration of the Ising model has a horizontal crossing if there exists 
a path made of edges between nearest neighbor sites going from the leftmost inner column 
to the rightmost one and visiting only plus spins. (Because of the relative position of the 
rectangle and the lattice, the vertical columns are made of sites in zigzag.) Let 7ih(V,H) 
be the probability at the critical temperature of such a crossing. One defines similarly 
a vertical crossing and the associated probability tt v (V,H). It is a well-known fact from 
percolation theory that, on the triangular lattice, a horizontal crossing on plus spins (or 
open sites) exist if and only if there is no vertical crossing on minus spins (closed sites). 
(One can convince oneself easily of this simple fact using a drawing.) Since plus and minus 
spins are equiprobable, this observation implies 717, (V, if) + tt v (V,H) = 1. In this paper we 
will ultimately be interested in the limit n h (r) = liniy^^^ r= ^ v /2H ^hiY, H) and similarly 
for tt v . The infinite lattice limit of the previous relation is known as the duality relation: 
7!7i(r) + 7r„(r) = 1. If iTh and n v are rotational invariants, the latter relation can be written 
as 7T/j(r) + 7!7j(l/r) = 1. This duality relation will hold for other (regular) lattices if the 
functions iTh and ir v are universal but its discrete equivalent (71^ (V, H) + ir v (V, H) = 1) hold 
strictly for neither the square nor the hexagonal lattices. 

The present paper discusses both a prediction for n h extending Cardy's approach and 
precise measurements of the function ir^ for 40 values of its parameter r. The agreement 
will be seen to be excellent, that is, perfect within statistical errors. The first section covers 
the theoretical prediction, the second the details of the simulation and the comparison of 
the data with the prediction. 
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2 A theoretical prediction based 
on conformal field theory 



2.1 Cardy's prediction for percolation 

Cardy's prediction for for two-dimensional percolation proceeds in two steps. He first 
identifies the probability ith with the difference of two partition functions with boundary 
conditions for the 1-state Potts model. He then uses the conformal field theory at c = to 
obtain an analytic expression for this difference. Here is a (very rapid) presentation of these 
two steps. 

The partition function Z(q) of the g-state Potts model on a finite rectangular domain 
is the sum over all configurations a of e~ /3i ^°") where H(a) = J^2(x,y)(^ ~ ^a(x),a(y))- The 
sum in H(a) runs over immediate neighbor pairs (x,y). This can be rewritten as Z(q) = 
Y1rP B {1 ~ p^ B ~ B ( R )q Nc ( R ) where p = 1 — e~^ J . The sum is over all subsets R of the set 
of edges of the lattice in the rectangular domain. The integer B counts the edges in the 
lattice, B(R) those in the subset R and N C (R) the clusters in R. If q — 1 (the value for 
percolation), this sum is 1 as desired. Let Z a p be the partition function of the Potts model for 
configurations whose spins on the left side of the rectangle are in the state a € {1,2,..., q}, 
those on the right in the state (3, and the others free. Cardy's first crucial observation is 
that ii h = {Z aa — Z a/3 )\ q=1 where a ^ (3. (The difference, done for a "generic" q, contains 
precisely the configurations that have a cluster intersecting the left and right sides.) The 
problem of calculating ir^ is therefore transformed into that of calculating partition functions. 

The possibility of calculating partition functions on finite domains with given boundary 
conditions also originates from works by Cardy (see for example |J). In the case of Z a p, for 
example, the four sides of the rectangle are submitted respectively to the boundary conditions 
a, free, (3 and free. Cardy argues that such a partition function is proportional to the 4- 
point correlation function ((j)(zi)(j)(z2)<p(zs)4>(z4)) in the conformal field theory associated to 
percolation whose central charge is c = 0. The Zi G C are the vertices of the rectangle in 
the complex plane and is the field that changes the boundary conditions in this theory. 
For percolation, this field is identified to 01,2 and it has conformal weight h — 0, a necessary 
condition for (0(zi)0(z2)0(z3)0(z4)) to be scale invariant. (The indices on 0^2 refer to 
the labels of Kac table. See 0.) The rules to find correlation functions in a conformal field 
theory are well-known and with this identification between partition functions with boundary 
conditions and 4-point correlation functions, the problem of calculating 7r^ amounts to solving 
an ordinary differential equation. 

Two obstacles appear in applying these ideas to the Ising model. First the Ising model 
is the q = 2 Potts model and the difference [Z aa — Z a p)\ q= 2 cannot be interpreted as the 
crossing probability since the factor q Nc ^ in the sum gives different weights to the various 
configurations that have a crossing from left to right. Second although the operator for 
the Ising model that changes the boundary state from free to a given state a is still 0^2, 
as in percolation, its conformal weight hip is now and its 4-point correlation is not 
anymore invariant under a conformal mapping z — > w but picks up the usual jacobian factors: 
(<j>{w 1 )<f>(w 2 )<p(w 3 )<j)(w i )) = (]\ i \w\z i )\~ hl - 2 ){(p(zi)(l)(z2)<p{z z )(t)(z i )). These prefactors (and 
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those additional coming from the presence of vertices along the boundary) seem to contradict 
the (strict) conformal invariance observed by simulation in ||. 

In || and |7j the probability tt^ of having simultaneous horizontal and vertical crossings 
was also obtained numerically Watts || was able to extend Cardy's argument to obtain 
a prediction that fits extremely well the data. His work is of particular interest to us 
as he is able to write 717^, again as a partition function but with more general boundary 
conditions. Whether a conformal boundary operator accomplishes the change between these 
more complex boundary conditions is not clear. Nonetheless the expression of 717^ as a 
partition function allows us to expect, Watts argues, that this probability is given by some 
4-point correlation function. He seeks it in the h = sector. 

2.2 The differential equation for tt^ 

Cardy's argument for percolation cannot be extended to the Ising model. One can still hope 
to relate crossing probabilities like i\h to 4-point correlation functions as Watts did for 717^ 
of percolation. If such a relationship exists, the choice can be narrowed to 4-point functions 
of the identity family as these are the only ones in the c = | conformal field theory that 
are invariant under conformal map z — > w like 7r^(r). An obvious objection will be that the 
4-point function of the primary field in the identity family is identically 1 (or a constant). 
For the calculation at hand it might well be that this primary field must be interpreted as 
one whose correlation functions satisfy only one of the differential equations corresponding 
to the two leading singular vectors. This milder requirement does not force the function to 
be a constant as will be seen immediately. 

The Verma module V/ C= i h=0 ^ of the Virasoro algebra has a maximal proper submodule 
M generated by two singular vectors, one at level 1, the other at level 6. The first of these 
is £-i|0) and the corresponding differential equation implies that the 4-point function is a 
constant. We shall drop this requirement. The other singular vector is 

(l^ - 10 L-\L-z + f L_x 2 L_ 2 2 - ^ L_ 2 3 + I L_! 3 L_ 3 

337 3185 2 381 2 1265 

6 144 6 2 18 

19309 T T 9005 T \ , . , . 

H L-iL- 5 L_ 6 c = ~, h = 0) . 

36 ° 12 7 1 2 ' 1 

If f(z) is the 4-point function with z={z\ — z^)[zz — Za)/{z\ — z^){z2 — 24), the differential 
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equation is 



^ (1 - 2 z) (686 z 2 (1 - zf + 73 z (1 - z) + 25) -^/(z) 

1 , , ^ „„„„ , ._.s rf 2 



+ — 2 (1 - z) (25141 z 2 (1 - zf - 2986 z (1 - z) - 171) — /(*) 



v _ „) 2 _ 

144 v ' v v ; v ' ' dz 2 

1 2 , 9 „ x _ d 3 



+— z 2 (I -zf (1-2 z) (208 - 3595 z (1 - z)) —f(z) 
27 dz A /-.\ 

1 r/ 4 ^ 

+-^ 3 (l-z) 3 (137- 737z(l-z)) — f(z) 

+ 10 (l-2z)z 4 (l-*) 4 ^/(z) 

^6 

+, 5 (1 - z f —f(z) = 0. 

This differential equation has three (regular) singular points at 0, 1 and oo. It is invariant 
under any permutation of these three points. (Invariance under z — > 1 — z is clear: invariance 
under z — > 1/z requires some work.) The exponents at any of these points are 0, | twice 
degenerate, |, | and |. The monodromy matrices around the three singular points are similar 
due to the symmetry of the equation but they cannot be diagonalized simultaneously. 

The cross-ratio z — (z\ — Z2){z% — z^)j{z\ — z 3 ) (z 2 — 24) is related to the aspect ratio 
r = y/3V/2H of the rectangle. If the four points zi, Z2, £3 and Z4 are chosen along the real 
axis at — |, —1, 1, |, then 
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1 + v^' 

A Schwarz-Christoffel transformation can be used to map the upper plane onto a rectangle 
with the images of the Z{ at the vertices. The aspect ratio is then given as 

_ K(l-k 2 ) 
T ~ 2K(k 2 ) 

where K is the complete elliptic integral of the first kind. The very short but wide rectangles 
(r = y/3V/2H — > + ) corresponds to z — > + , the tall and narrow (r — > +00) to 2 — > Y~ and 
the square to r = 1, 2 = |. The function r(z) has the property r(z) = r(~) and the symmetry 
2; — ► - of the differential equation is thus welcome. If 717, does not depend on the relative 
angle between the rectangle and the lattice, that is if Hh is a rotational invariant, then the 
duality relation implies 71^(1) = | and it can be put in the form (| — 7Tft(r)) = — (| — 7Tfc(-)). 
Fortunately the function r(z) is such that r(z) = l/r(l — z) and the duality simply states 
that the function f(z) = | — 7T/ l (r(z)) is odd with respect to the axis 2 = |. An odd subspace 
of the solution space of the differential equation (odd with respect with z — ~) exists due to 
the symmetry z — > 1 — z and it is of dimension 3. 

We explored several paths to cast the solutions of eq. ([!]) into analytically tractable forms. 
One of them was to write the lhs of eq. ([!]) as rii<j<6( z (-'- ~~ z )f l tL ^ e Watts did. But there 
is no real solutions for the Oj's in the present case. The most natural path however is the 
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screening operator method (Q, see also ||). The pertinent field is 2 ,3 (with Kac's labels) 
and the 4-point correlation calls for three contour integrals (a charge Q+Q 2 _ must be added 
at infinity to assure neutrality). Integral representations of six linearly independent solutions 
can be obtained in this straightforward (but probably tedious) way. The main problem is 
therefore whether there is sufficient physical information on / to fix the linear combination. 
As argued above the odd parity of / reduces the space of solutions to a three-dimensional 
subspace. The condition /(0) = — \ (i.e. 7^(0) = 1 — 71^(1) = 0) is one further linear 
constraint. The function f(z) is monotone increasing but this condition will restrict / to 
an open set of the 3-d subspace rather than decrease the dimension. The numerical data 
presented below indicate that tt^z) — > z& as z — > + . This is rather striking in view of the 
two- fold degeneracy of the exponent ~. Two solutions associated to this exponent can be 
chosen to behave as z^ and z^ log z for small z and the latter, if present in /, should dominate 
the former whenever z is close to 0. The fact that it is not seen in the simulation could be 
interpreted physically as a manifestation of the power law behavior of critical correlation 
functions at short distance. Imposing that the behavior in zg logz be absent of / would add 
one linear constraint. With all these constraints we would still be left with a one-dimensional 
subspace in the space of odd solutions. We have not found any further constraints to fix 
completely the function /. This is why we resorted to a numerical fit (see Paragraph |3.2| ). 

As we were trying to solve analytically the differential equation, Marc-Andre Lewis sug- 
gested to us to look for a solution of the form ez d 2Fi(a, b, c; z), with a, b, c, d and e constants, 
that would be odd around z = \ and reproduce the asymptotic behavior of the data. Cardy's 
prediction for percolation is of this form and the suggestion is natural in this sense. Such a 
function exists but it does not satisfy the differential equation. However it follows so closely 
the data of || (the worst gap is 1%) that we decided to improve these measurements to an- 
swer the question: which of the hypergeometric function or of the solution of the differential 
equation, if any, describes the data. 

3 Improved measurements of 

3.1 Finite size effects and power law behavior 

The probabilities 7ih, ti v and tt^ at 81 values of r (e [0.136, 7.351]) were measured in || for the 
three regular lattices on rectangles containing around 40000 sites. For these sizes, departure 
from the duality relation is small but still noticeable for the square and the hexagonal lattices. 
No attempt was made there to use various sizes in order to reduce finite-size effects. For 
the new runs to be presented here we chose to concentrate on the triangular lattice and use 
various sizes to approximate the function 7th (and the other two, n v and irhv) in the limit 
when the mesh goes to zero. 

A power law for the finite-size behavior of critical data is an accepted hypothesis and our 
measurements rest upon it. It states that, for sufficiently large size, 

\Tr h (V,H)-Tr h (r)\naVP 

with (3 a negative constant and r = \/3V/2H. We shall use several linear sizes of the form 
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V = 2 l V and H = 2 l H . Writing n h (i) for ^ h {2 l V G , 2 l H ) and supposing that ir h (i) is 
decreasing, we can write 7r h (i) - 7T h (i + 1) « aV^2^(l - 2^) and therefore express ^(r) as 
7r^(z) — oVq 2 l P. To determine the constants a and /? requires at least 3 rectangle sizes as 
only the differences (^(z) — ith(i + 1)) can be used. 



r 


Vo 


Ho 


tt*(0) 




7T ft (2) 


^(3) 


*A (4) 


0.1443 


4 


24 


0.029265|21 


0.025583|20 


0.023806|19 


0.022949|19 


0.022538|18 


0.9897 


16 


14 


0.502797|46 


0.499753|63 


0.498280|63 


0.497645|63 


0.497321|60 


6.928 


32 


4 


0.978259|17 


0.977944| 15 


0.977865|13 


0.977834|16 


0.977859| 18 



Table 1: ir/j for three aspect ratios r and five sizes. 



What are the right rectangle sizes and what is the required precision on each 7Tft,(^) ? For 
the two extreme rectangles that we are planning to measure, r wi d e = = 2^a ~ 0.1443 

and r ta ii = 2 T ~ 6.928, and for a rectangle close to a square r sq = ~ 0.9897, we 

obtained 7Tft(i), z = 0, 1,2,3,4. Each increment corresponds to an increase by a factor of 2 
of the linear size. For example 71^(0) was measured on a rectangle of 4 x 24 sites and 71^(4) 
on 64 x 384 sites for the wide rectangle. The results appear in Table [I]. The samples were 
large, at least 250 x 10 6 configurations. The digits after the vertical bar gives the statistical 
error on the digits just before; for example, the first element in the table (0.029269|21) 
means that Tf h (0) is 0.029269 with the 95%-confidence interval being [0.029248,0.029290]. 
The differences between 7r(3) and 7r(4) are however small. In fact the monotonicity of 7r^(i) 
for rtau = 6.928 is broken for i — 4 even though the error bars do allow for the power law 
to hold. Larger samples would definitely be required for the large lattices. Fortunately the 
precision on the measurements for i = 0,1,2 and the fact that the power law seems to hold 
for very small sizes (4 sites in one direction!) allow for good estimates of 7r^(r) without 
these larger lattices. Table |2| shows estimates of n h for r wide , r ta n and r sq using the power 
law hypothesis and a subset of the measurements of Table |l|. The notation i—j means that 
7T/i(i), 7th(i + 1), • • • , fth{j) were used to obtain ^(r). Using only the three smallest lattices, 
the three largest or the five ones lead to estimates 7r^(r) that differ by less than 4 units on the 
fourth significant digits. We therefore decided to use only three sizes for each r considered 
and choose the pairs (Vo, H ) in such a way that V ,H > 4 and that V H > 96. All samples 
were larger or equal to 10 8 . Tables |3] gives the results for n h , n v and ir hv at 40 values of 
r e [0.1443,6.928]. 



r 


0-2 


0-3 


0-4 


1-3 


1-4 


2-4 


0.1443 


0.02215 


0.02215 


0.02216 


0.02215 


0.02216 


0.02216 


0.9897 


0.4969 


0.4971 


0.4970 


0.4972 


0.4970 


0.4970 


6.928 


0.9778 


0.9778 


0.9778 


0.9778 


0.9778 


0.9777 



Table 2: Estimates 7r h using subsets of available data. 



There is always, on a finite lattice, the problem of determining r from the numbers V and 
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H. The two simplest choices are the aspect ratios of the smallest or the largest rectangles 
that include the sites considered and only those. Even though it is not a natural choice, 
the aspect ratio of the tallest and narrowest rectangle would be another convention. The 
method we have used to determine 717, (r) overcomes this imprecision due to convention. For 
any convention the aspect ratio for a rectangular subset of the triangular lattice will be 
^(V + A v ) / (H + A H ) with Ay and A H dependent on the convention but independent of V 

and H. The ratio r at which 717, is measured is therefore lim^oo ^(2 l V + Ay) / (2 l H + A H ) = 

^Vq/Hq, independent of Ay and A H , that is independent of the convention. This is another 
advantage of using several lattices for a given r. 

Even though 717, (and tt v ) is invariant under rotation, finite-size effects are not. The 
differences between fc h (0) = 0.02927, w h (l) = 0.02558 and f} h (2) = 0.02381 for r = 0.1443 
are much bigger than those between fc v (0) = 0.02174, fc v (l) = 0.02205 and fr v (2) = 0.02214 
for 1/r = 1/6.928 = 0.1443. The statistical errors on 71^(0.1443) and 7r„(6.928) are therefore 
different even if both numbers turn out to be very close (0.02215 and 0.02218). In the worst 
cases the 95%-confidence interval amounts to less than 2 units on the third significant digit 
(e.g. 7^(0.1443) = 0.02215 ± 0.00015). At the center of the range of r the error on 7Tft(r) 
decreases to 4 units on the fourth digit and it is even smaller for large r. Over the whole 
range it is smaller than 4 x 10 -4 for both and 7r„. 

The improvement upon previous measurements found in || can be checked easily. Among 
the 40 values of r used there are nine pairs ((V a , H a ), (V&, H b )) such that 

Ta 2 H a \2 H b ) " r b 

The pairs (r a , l/r ) are those corresponding to the following lines of Table |3|: (1, 40), (5, 36), 
(7, 34), (8, 33), (10, 31), (11, 30), (12, 29), (18, 23), (20, 21). The measurements for these pairs 
should satisfy f\ h {r a ) = 7t v (r b ), f\ v {r a ) = n h (rb) and f\ hv {r a ) = f} hv {r h ) within statistical errors. 
This turns out to be the case. (See Figure |1|.) There are in total 27 independent comparisons. 
Their relative errors is always less than 2 x 10~ 3 . The largest occur when the quantities f} 
being compared are themselves very small, like 7T/i(-^|) = 0.02215 and ^(^8) = 0.02218. 
In all cases the absolute value of these differences are less than 5 x 10~ 4 . Not only are these 
variations small, they are of both signs. This fact is a further indication that the power 
law hypothesis provides a very good approximation. Suppose indeed that, at the sizes used, 
correction terms are required: \^h{i) ~ 7r ^( r )l ~ aV^il + y + . . . ). These new terms would 
lead to a systematic error that is not seen here. 

3.2 A prediction for 7Th for the critical Ising model 

The first easy comparison between the theory developed in the first Section and the new 
data lies in the asymptotic behavior of 7r^ as r approaches and +oo. If ir^ is a solution of 
the differential equation (|1|), then \ogn h (r) — > —Xn/r as r — » and log(l — 7f h (r)) — > —Xirr 
as r — > +oo with A one of the exponents. Using the ten extreme values of r on each side of 
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Figure 1: Differences n h (r a ) - n v (r b ) and n hv (r a ) - n hv (r b ) for several pairs (r a , r b = r a l ). 

the measured interval, we obtain for 7r/j 

log 7r h (r) — > -0.166487T- 

i — >o r 

log(l — TTh(r)) — > — 0.166577rr. 



The slopes obtained using the data for ir v are 0.166477T and 0.166547T. These numbers are 
very close to the exponent |, the smallest non- vanishing exponent of the differential equation 
(0). This is remarkable! Clearly A is to be interpreted as a critical exponent of the Ising 
model. But none of the usual exponents of the Ising model contains the prime number 3 in 
their denominators and scaling laws involve only products and integral linear combinations 
of these exponents. If this new critical exponent can be deduced from the usual ones, it will 
not be by traditional scaling laws. 

The second test is to obtain a solution of the differential equation ([!]) that describes the 
data. Let fi, i = 0, 1, . . . , 5, be a basis of solutions for equation (jl]) defined by their behavior 
at z — |: fj^ (z = |) = 8ij. Since f(z) — \ — ir h (r(z)) is odd with respect to z — |, it lies 
in the subspace of functions of the form af\ + bf 3 + cf 5 . We determine the constants a, b 
and c by requiring that L = — ^h( r i) — /(- 2 ( r i))) 2 ) 5 be minimum. The sum is over 

the 40 data. The three solutions f\, fs and f$ were obtained numerically. Both Matlab 
and Mathematica give similar fits. These softwares have internal parameters controlling the 
required accuracy of the integration. These parameters can be pushed to a point where 
stronger requirement does not lead to any significant improvement on the minimum of L. 
Figure ^ has been drawn using the values a, b and c obtained for control parameters beyond 
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Figure 2: The prediction log7T/i/(l — 7i~h) as a function of log r together with the 40 measure- 
ments TTh(ri). 

this point. The largest among the differences ||— Tth{ r i)— f{ z { r i))\> i = 1, ■ ■ ■ , 40, is 3.6x 10~ 4 , 
smaller than the statistical error, and the standard deviation is 1.5 x 10~ 4 . Similar results are 
obtained for tt v . The agreement is therefore excellent. As a comparison it is instructive to 
redo Cardy's calculation testing his prediction for percolation. Since the publication of [l|, 
better data were obtained for percolation by sites on a square lattice for 81 rectangles with at 
least 10 6 sites [0. The samples contained over 10 6 configurations. For these, the statistical 
errors are approximately 10~ 4 at the extremities of the interval (r G [0.142, 7.067]) and 10~ 3 
in the middle. (At the extremities these data are therefore more precise than the present 
ones for the Ising model and at the center they are less precise.) The largest departure 
from duality is 4 x 10~ 4 , almost exactly what is seen in Figure [I] for the present data. 
Comparing his prediction 7r^ orco (r(z)) = 3T(^)zs 2 F 1 (^, |, |, z)/T(|) 2 with the data leads to 
maxims! \TTh Crc °{r(zi)) — Tih{ r d\ = 7-8 x 10~ 4 and to a standard deviation of 4.2 x 10 -5 . 
These results are similar to those just reported. 

We mentioned earlier the possibility of describing the data with a hypergeometric func- 
tion. The function 

is odd around z — |, when shifted by |, and behaves as and (1 — 2) a when z — > + and 
z — > 1~. Several data are now more than 2 x 10 -3 apart from the corresponding values of 
g, a gap barely visible on a figure, but clearly out of any reasonable confidence interval. We 
may accept the solution of the differential equation as an analytic prediction for 7th but we 
must reject the function g(z). 
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Measurements were also made of the probability tt^. As its behavior at z = and z = 1 
is also described by the exponent h, one may hope that the even subspace (around z = |) 
of the differential equation may contain a solution matching the data. This is not the case. 
The best fit in this subspace lies up to 6 x 1CT 3 away from the data, an unacceptable gap. 
(This disagreement has an advantage. It shows that the success of the fit for 7T/i is not a 
consequence of the large freedom that a three-parameter fit gives.) Watts used the third 
singular vector to describe %hv for percolation. In the Verma module V c=0t h=o, this vector is 
at level five, leading to a differential equation of order 5. However the third one in V c= i h=Q 

is at level 11 and the associated operator P(z, 4-) is a polynomial of order 11 in j-. It can 

be cast as 

p ^)= £ ((z(i-zr ft ( 2 (i-z))^ 

l<i<ll, i odd 

+ J2 (l-2z)((z(l-z)y- 1 p i (z(l-z))^- 

2<j<10, i even 

with 

Vl {u) = -±u (-484 - 20465 m + 120702m 2 + 134456 m 3 + 326536 m 4 ) 
p 2 {u) = -± (-484 - 65267 m + 555942 m 2 + 2442422 u 3 + 9038782m 4 ) 
p 3 (u) = (44318 - 1002950 u - 5132553 u 2 - 12317152 u 3 + 245463307 u 4 ) 
p A {u) = ^ (-412839 - 8111249 u - 42237350 u 2 + 749236363 m 3 ) 
p 5 ( u ) = - (3492203 + 28198986 u - 1347384726 u 2 + 5106251212 m 3 ) 
p 6 (u) = — Yg^4 (368339 - 27634444 u + 153070553 u 2 ) 
p 7 ( u ) = (554551 - 9066926 m + 29200567m 2 ) 
p 8 (u) = g (-11842 + 75835 u) 
p 9 (u) = -11 (-757 + 3457 m) 

Pio(«) = 
Pu(m) = 1. 

The differential equation P(z, j-)f(z) = is again symmetric under any permutation of 
the three regular singular points 0, 1 and oo. The exponents (with their degeneracies) are 
0(2), |(2), |(1), 1(1), |(2), |(1), an d 6(1). The even subspace of the differential equation 
is of dimension 6 and the fit will therefore contain 6 parameters. Using the data of the last 
column of Table |3] and numerical integration of the differential equation, we obtain a best 
fit (in the same sense as the one used for tt^) that has a largest difference of 2.7 x 10 -4 and 
a standard deviation of 1.3 x 10~ 4 . These are excellent results well within the experimental 
windows. However we have tried to fit similarly the function h(z) = kz&(1 — z)e with the 
constant k chosen such that h(z = |) = 7Ch v (z = |). While this function is not a solution, the 
6-dimensional subspace of (numerical) even solutions contains a solution h that approaches 
it extremely well (namely ma.x z& [o^ \h(z) — h(z)\ ~ 10~ 6 ). Consequently the large dimension 
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of the even subspace allows for functions that are not solution to be fitted within statistical 
errors and the fit for ir^ is much less convincing than the one above for tt^ or than Watts' 
prediction. 
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Table 3: The measurements fc h , fr v and 717^ 
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